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Abstract 

We apply the perturbative QCD factorization theorem developed recently 
for nonleptonic heavy meson decays to the radiative decay B — * K*"f. In this 
formalism the evolution of the Wilson coefficients from the W boson mass 
down to the characteristic scale of the decay process is governed by the effec- 
tive weak Hamiltonian. The evolution from the characteristic scale to a lower 
hadronic scale is formulated by the Sudakov resummation. Besides comput- 
ing the dominant contribution arising from the magnetic-penguin operator 
O7, we also calculate the contributions of four-quark operators. By fitting 
our prediction for the branching ratio of the B — * K*"f decay to the CLEO 
data, we determine the B meson wave function, that possesses a sharp peak 
at a low momentum fraction. 
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I. INTRODUCTION 



The observation of the decay B — > K*j five years ago [Q] opened a new era for particle 
physics, since the penguin structure of the electroweak theory was probed for the first time. 
Soon after this observation, the inclusive B — > X s 7 decay was also established. The 
updated branching ratios for both decays are (4.2 ± 0.8 ± 0.6) x 10 -5 and (3.14 ±0.48) x 
10 -4 0J|, respectively. In recent literatures, the inclusive decay B — > X s ^ receives more 
attentions than the exclusive mode B — > K*^. The branching ratio and the photon energy 
spectrum of the B — > X s j decay have been used to constrain the parameter-space of the new 
physics beyond the standard model. The preference of studying the inclusive mode is clear 
in that its hadronic dynamics is much easier to handle as compared to the exclusive mode. It 
has been shown that, to the leading order in I /Mb, Mt> being the b quark mass, the branching 
ratio B(B — > X s j) is given by the branching ratio B(b — > 57) of the corresponding quark-level 
process. Furthermore, the sub-leading l/M^-corrections can be parametrized systematically 
using the heavy quark effective theory ||. 

The hadronic dynamics of the exclusive decay B — > K*^ is much more complicated. 
Specifically, one has to deal with the soft dynamics involved in the B and K* mesons. 
Since the final states are light compared to the decaying B meson, one anticipates that 
perturbative QCD (PQCD) is applicable to this process because of the large energy release. 
In fact, the PQCD formalism based upon factorization theorems, which incorporates the 
Sudakov resummation of soft-gluon effects, has been developed for sometime |7|-£|, and 
applied to semileptonic f/]||, inclusive-radiative || and nonleptonic B meson decays |10 



This approach is so far rather successful. In this article, we shall extend this formalism to 
penguin-induced exclusive processes, such as the B — > K*^y decay. The satisfactory result 
with repect to this complicated process, as presented later, provides further confidence on 
the validity of the PQCD approach to B meson decays. 

According to the PQCD factorization theorem, the branching ratio of a heavy meson 
decay is expressed as the convolution of a hard subamplitude with meson wave functions. 
The former, with at least one hard gluon attaching to the spectator quark, is calculable 
in the usual perturbation theory, while the latter must be extracted from the experimental 
data or derived by nonperturbative methods, such as QCD sum rules. Since the K* meson 
wave function has been known from the sum rule analyses [ |TT| , |T2"|| , the B —>■ K*j decay 
is an ideal process from which the unknown B meson wave function can be determined. 
With the B meson wave function obtained here, we are able to make predictions for other 
B meson decay modes, especially for the charmless decays. We shall show that by fitting 
our prediction for the branching ratio B(B — > K*^) to the CLEO data, a B meson wave 
function with a sharp peak at the low momentum fraction is obtained. 

Comparing our work to others, we remark that most of the previous studies on this decay 
focus only on the contribution of the magnetic-penguin operator O7, and their approaches are 
based upon quark models. To ensure that other contributions are indeed negligible, we also 
calculate the contribution by the current-current operator O2 , which arises through the four- 
point b — > sg*7 coupling with the off-shell gluon reabsorbed by a spectator quark. Although 
such a contribution has been calculated before |R|, it is however obtained in a naive PQCD 



framework which does not include the Sudakov resummation and the renormalization-group 
(RG) analysis |y|. As a result, the predictions of Ref. 0] are sensitive to the choice of 
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the renormalization scale \i. In our more completed approach, the sensitivity to /x can be 
avoided. We find that, due to certain cancellations, the contribution by 2 turns out to be 
rather small. 

We also like to comment on a different viewpoint based on the overlap integral of meson 



wave functions 14], which showed that the diagrams without any hard gluon dominate over 



those we will be considering here. We shall argue that the observation in |14] is due to an 
underestimation on the value of the strong coupling constant and a choice of a flat B meson 
wave function. If evaluating the coupling constant at the characteristic momentum flow 
involved in the decay process and employing a sharper B meson wave function, these 
higher-order contributions may become comparable to the leading-order ones. Hence, the 



approach in JL4J] does not seem to be self-consistent. In our approach, the diagrams without 
hard gluons do not contribute under our parametrization of parton momenta. Therefore, 
the diagrams with a hard gluon attaching to a spectator are leading in our analysis. 

This article is organized as follows: In Sect. II, we write down the effective Hamiltonian 
for the B — > K*^ decay. In particular, the current- current operator O2 and the magnetic- 
penguin operators O7 are identified as major sources of the contribution. We then calculate 
the 02-induced b — ► sg*^ vertex, keeping the gluon line off-shell. In Sect. Ill, we derive the 
factorization formulas for the B —>■ K*^ decay, which include contributions from the various 
operators. The numerical result is presented in Sect. IV, where the contribution of each 
operator is compared. Sect. V is the conclusion. 



II. EFFECTIVE HAMILTONIAN 

The effective Hamiltonian for the flavor-changing b — > s transition is now standard, which 
is given by p"5| . p"6 



with 



H eS (b - s 7 ) = -^V*EWW C 1 ) 



X = {SiCj)v-A{Cjbi) V - A 

2 = (SiCi) V -A(Cjbj)v-A 

3 = (sibi) v _ A J2(qjqj)v~A 

1 

4 = (sibj) V - A J2{QjQi)v-A 

<? 

5 = {Sibi)v-A^2{mj)v+A 

<? 

Or = £^SiCT^(m s P L + m b P R )biF^ 

8 = ^s l( r^(m s P L + m b P R )T^G a ^ (2) 
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i, j being the color indices. In the PQCD picture, the lowest-order diagrams for the B — > 
K*^ decay arising from operators 2 , O7 and 8 are depicted in Figs. 1-3. It is not hard 
to see that contributions other than the above are negligible. For example, contributions 
depicted in Fig. 4 are very suppressed, although they are of the same order, 0(eGpa s ) |17| , 
as those of Figs. 1-3. We draw this conclusion from a previous experience with B — > D*'-/. 
Indeed, from a diagram similar to Fig. 4 (with s replaced by c, and the q on both sides 
replaced by d and u, respectively), one has obtained B(B — ► D*j) = 10~ 6 Since 
fx* ~ Id*, an d VtbVts is comparable to V c bV u d, it is clear that the diagram in Fig. 4 
gives B(B — > « 10 -6 x Cf, 4)5)6 s « 10~ 9 — 10~ 10 , and is thus negligible. There is 

still one more type of contributions of the same order as shown in Fig. 5. By an explicit 
calculation, one can show that such contributions merely give corrections to the Wilson 
coefficients C3-CQ occuring in Fig. 4. Since the matrix elements of Fig. 4 and Fig. 5 have 
identical tensor structures, the branching ratio contributed by the latter figure behaves like 
10~ 6 x Cf x (^) 2 ~ 1CT 9 , which is also negligible. Finally, in Sect. IV, we shall see that 
the contribution by Og is negligible as well. 

Before implementing the PQCD formalism to evaluate the contributions of Figs. 1-3, we 
compute the four-point b — > sg*^ vertex. Our calculation essentially generalizes the work 



by Liu and Yao JOJ to the off-shell gluon case [2D]. We first perform a Fierz transformation 
on 2 , i.e., 

2 = (SiCi) V ^ A (Cjbj)v-A = {Sibj)v-A(CjCi) V -A. (3) 

The b(p) — > s(p')j(ki)g*(k2) vertex is expressed as 

I au = C 2 V t :V tb u(p')^Y(l - 75 )T a u(p) W M (^i)> ( 4 ) 
with the structure-tensor 

Ipvp = A\€ pl j pcr ki + A 2 e piUpa k 2 + -A^ppaT k\ k 2 k\ v 

+ A^e vpuT k x tpparki k 2 k 2v ~\~ A§£ vpaT k^ k 2 k\fi- (5) 

Clearly, the form factor Aq can be discarded because of k\ ■ e{k\) = 0. From the requirement 
of gauge invariance, i.e., k^I^p = and k 2 I pup = 0, we have 

A 2 + A 4 h -k 2 = 

A x + A 3 k! ■ k 2 + A 5 k 2 2 = 0. (6) 

The invariance of l pvp under the interchanges k\ <-> k 2 and \i <-> v further require A 3 = —A 4 . 
With the above relations, I^ up is simplified into 

Ifwp = ^4 [ki • k 2 e pvpa {ki — k 2 ) + e I/paT k 1 k 2 k 2p — e jip(TT k 1 k 2 k\ u ] 

+ A<z [tp P (jTkik 2 k 2l , k 2 t pvpcj k-^j , (7) 



with 



2\/2 

^4 = ^e^G F / n (M 6 2 ), (8) 
A 5 = -^eg s G F \l w (M 2 )-I 20 (M 2 )} , (9) 
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and 



Iab(m 2 ) = I dx I dy— X V — r-— . 

Jo Jo x{± — x)k 2 + zxyki ■ k 2 — m z + is 

Carrying out the Feynman-parameter integrations, we obtain 
/ii(m 2 ) = 



jgrations, we obtain 

1 r l , m 2 — x(l — x)k 2 
+ dx- 



2Q 



xQ 4 



-In 



m 2 - x(l - x)(k% + Q 2 



ho(m 2 ) - I 20 (m 2 ) = 



m 2 - ar(l - x)(A; 2 + Q 2 ) 



m 2 — x(l — x)k\ 



(ix 1 In 



m 2 — x(l — x)k\ 



171 



&(Q 2 + k 2 2 -4m 2 )J^ 



(10) 



(11) 



(12) 



with Q 2 = 2h 
Iu(m 2 ) and hoim 



k 2 , and (3 = 

2\ T (^,1 



l-4m 2 /{k% + Q 2 



— hoim 2 ) have smooth limits as m 2 



As a side remark, we note that both 
0. Hence, we can safely neglect 
sg*7 vertex, since their Wilson coefficients, C4 



the contributions of 4 and 6 to the b 
and C*6, are much smaller than C 2 - 

Having determined the b — > sg*^ vertex, we attach the off-shell gluon line to the spectator 
quark to form the B — > K*^ amplitude as shown in Fig. 1. As mentioned, this amplitude is 
of the order ea s GV, the same as the amplitudes induced by O7 and 0$ depicted in Figs. 2 
and 3, respectively. In the next section, we shall compute the contribution of each operator 
to the B — > K*^f decay using the PQCD factorization theorem. 



III. FACTORIZATION FORMULAS 

In this section we first review the PQCD factorization theorem developed for nonleptonic 
heavy meson decays |H| , and then extend it to the B — > K*j decay. Nonleptonic heavy 
meson decays involve three scales: the W boson mass Mw, at which the matching conditions 
of the effective Hamiltonian to the original Hamiltonian are defined, the typical scale t of a 
hard subamplitude, which reflects the dynamics of heavy meson decays, and the factorization 
scale 1/6, with b the conjugate variable of parton transverse momenta. The dynamics below 
1/b is regarded as being completely nonperturbative, and parametrized into a meson wave 
fucntion (j)(x), x being the momentum fraction. Above the scale 1/b, PQCD is reliable 
and radiative corrections produce two types of large logarithms: \n(M w /t) and ln(tfe). The 
former are summed by RG equations to give the evolution from My/ down to t described 
by the Wilson coefficients c(t) . While the latter are summed to give the evolution from t to 
l/b- 

There exist also double logarithms ln 2 (P6) from the overlap of collinear and soft di- 
vergences, P being the dominant light-cone component of a meson momentum. The re- 
summation of these double logarithms leads to a Sudakov form factor exp[— s(P, b)], which 
suppresses the long-distance contributions in the large b region, and vanishes as b = 1/A, 
A = Aqcd being the QCD scale. This factor improves the applicability of PQCD around 
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the energy scale of few GeV. The b quark mass scale is located in the range of applicability. 
This is the motivation we develop the PQCD formalism for heavy hadron decays. For the 
detailed derivation of the relevant Sudakov form factors, please refer to [|7||§. 

With all the large logarithms organized, the remaining finite contributions are absorbed 
into a hard b quark decay subamplitude H(t). Because of Sudakov suppression, the per- 
turbative expansion of H in the coupling constant a s makes sense. Therefore, a three-scale 
factorization formula is given by the typical expression, 



c(t) (8 H(t) ® <p(x) <g> exp 



-s{P,b)-2f ^ q (a s (p)) 
Ji/b fi 



(13) 



where the exponential containing the quark anomalous dimension 7 g = —a s /7c describes the 
evolution from t to 1/6 mentioned above. The explicit expression of the exponent s is can be 
found in [ID]. Since logarithmic corrections have been summed by RG equations, the above 



factorization formula does not depend on the renormalization scale \i |T0j. Our formalism 
then avoids the sensitivity to fi that appears in [13" . 



We now apply the three-scale factorization theorem to the radiative decay B — > K*j, 
whose effective Hamiltonian has been given in the previous section. As stated before, only the 
operators O2, O7 and Og are crucial, to which the corresponding diagrams and hard subam- 
plitudes are shown in Figs. 1-3 and Table I, respectively. We write the momenta of B and K* 
mesons in light-cone coodinates as Pb = (Mb/v / 2)(1, 1,0t) and P% = (Ms/\/2)(l, r 2 , Ox), 
respectively, with r = Mr* /Mb- The B meson is at rest with the above choice of momenta. 
We further parametrize the momenta of the light valence quarks in the B and K* mesons 
as ks and kx, respectively, ks has a minus component k B , giving the momentum fraction 
xb = kg/P^, and small transverse components ksr- kx has a large plus component k K , 
giving %k = k^/P^, and small k KT . The photon momentum is then P 7 = P B — P K , whose 
nonvanishing component is only P~. 

The B — > K*^j decay amplitude can be decomposed as 

M — e*- e* K ,M s + ie, p+ ^e*£,M p , (14) 

with e 7 and ex* the polarization vectors of the photon and of the K* meson, respectively. 
Note that we have neglected the structure (P 7 • e* K *)(PK • e*)/(P 7 • Pk) which should come 
together with e* • e* K ,. This is due to our choice of the frame which gives Pk • e* = 0. From 
Eq. (|14D, it is obvious that only the K* mesons with transverse polarizations are produced 
in the decay. 

The total rate of the B — > K*^ decay is given by 



2 

'' (|M 5 | 2 + |M P | 2 ) . (15) 



8ttM b 

We can further decompose M s and M p as 

M i = M^ + Mij + Mi, (16) 

where % = S or P, and the terms on the right-hand side represent contributions from 
operators 2 , O7, and 8 , respectively. 
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In the following, we write the factorization formulas for M\ in terms of the overall factor 

T (0) = ^V*V tb C F M B . (17) 



The amplitudes contributed by 2 is written as 

Mf = T^ '- dx dy dxsdxx \ bdb(j) B (x B )(f)K*(xK) 
3 Jo Jo Jo Jo 

xa s (t 2 )c 2 (t 2 ) exp[-S(x B , x K , t 2 , b, b)] 

x [(1 — r 2 + 2txk + 2xb)v ~~ ( rx i< + 3xb)(1 — x)] 

(1 — r)(l — r 2 )!^! 
X xy(l-r 2 )x x M|-M c 2 



H 2 (Ab,J\Bl\b) 



(18) 



with 



1 ,1/A 

dxsdxx / bdb<pB(xB)4>K*{ x K) 
jo 



f dx j X ' X dy 
6 Jo Jo 

xa s (t 2 )c 2 (t 2 ) exp[-5'(a; B , x K , £ 2) &, &)] 
x "((l-r)(l 

- (r(l + r)x^ + (3 - r)x B ) (1 

(l-r 2 )x^x rr/ ,, /TTT^v 



r 2 ) + 2r 2 xx + 2xb) V 



x)\ 



(19) 



A 2 



x K x B M 2 B , 
x K x B Ml - 



1 



x 



1 - r 2 )x„M 2 + 



M 2 



x 1 



x) 



max(A, 



Bl\,l/b). 



(20) 



To arrive at the above expressions, we have employed Eq. (|T^) for the charm loop integral, 
instead of Eqs. ( |TT] ) and ([12]). The variable bs [bx), conjugate to the parton transverse 
momentum k B T {kxT}-, represents the transverse extent of the B (K*) meson; t 2 is the 
characteristic scale of the hard subamplitude 



H 2 (Ab, J\Bl\b) = Ko(Ab) - K (J\B$\b) , B\ > 



K (Ab)-i^-H^( 



Bl\b) 



B z 2 <0 



(21) 



which comes from the Fourier transform of the corresponding expressions in Table I to 
the b space. Note that, for simplicity in the notation, we have used H 2 to denote hard 
subamplitudes in both M 2 and M 2 '. In Table I, these two amplitudes are distinguished. 

There are two diagrams, Figs. 2(a) and 2(b), associated with the operator O7, where the 
hard gluon connects both quarks of the B meson or those of the K* meson. The amplitudes 
from the two diagrams are 



7 



M 7 S = -M 7 P = r<°>2 / dx B dx K / b B db B b K db K <j) B (xB)<pK*(x K )(l - r 2 ) 
Jo Jo 



X 



rHi a \Ab K , y/\Bf\b B , y/\Bf\b K )F 7 (t 7a ) 



+ [1 + r + (1 - 2r)x K ]H?\Ab B , C 7 b B , C 7 b K )F 7 (t 7b ) 



B 2 = (x B - r 2 )M 2 B , C 2 = x K M 2 B 



(22) 



with 



t 7a = max(A, yJ\B%\, l/b B , l/b K ) , 
t 7b = max(A, C 7 , l/b B , l/b K ) , 

where the function F 7 denotes the product 

F 7 (t) = a s (t)c 7 (t) exp[-S(x B ,x K ,t,b B ,b K )] . 

The hard functions 

Hi a) (Ab K , ^/\Bj\b B , \[\Bj\b K ) = 

B 2 > 



(23) 



(24) 



K (Ab K )h(J\B 2 \b B ,J\B 2 \b K 



K (Ab K )h'(J\B 2 \b B , J\B 2 \b K ) , B 2 < 



(25) 



with 



h = 6{b B - b K )K (^\B 2 \b B )I (^\B 2 \b K ) + (6 B 



and 



6{b B - b K )H { G 1 \^\E 2 \b B )U^\Bj\b K ) + (6 B <- 6a:) 



H^\Ab B ,C-rb B ,C 7 b K ) = K (Ab B )h(C 7 b B ,C 7 b K ) , 



(26) 



(27) 



are derived from Figs. 2(a) and 2(b), respectively. The relation Mf = —Mf reflects the 
equality of the parity-conserving and parity-violating contributions induced by 7 . 

Four diagrams, Figs. 3(a)-(d), are associated with the operator 8 , where the photon is 
radiated by each quark in the B or K* mesons. The corresponding amplitudes are 



Mg = -r (0) - / dx B dx K b B db B b K db K (f) B (x B )(f) K *(x K ) 

O JO JO 

x |(1 - r 2 + x B )(rx K + x B )H { 8 a) (Ab K , B 8 b B , B 8 b K )F 8 (t 8a ) 

+ [(2 - 3r)x K - x B + r(l - x K )(rx K - 2rx B + 3x B )] 

xH^ b \Ab B ,C 8 b B ,C 8 b K )F 8 (t 8b ) 
+ (1 + r)(l - r 2 )[(l + r)x B - rx K ] 

xHi c) (J\A^\b K , D 8 b B , D 8 b K )F 8 (t 8c ) 



,1/A 



8 



with 

A' 2 

r 2 

t&a 
hb 

t 8 c 



-[(1 - r 2 )((l - r 2 )(2 - x K ) + (1 + Sr)(2x K - x B )) 
+2r 2 x K (x K - x B )] 



x (1 - r 2 + x B )(rx K + x B )W\Ab K , B 8 b B , B 8 b K )F 8 (t 8a ) 



+ [(2 - 3r)x K - x B - r(l - x K ){rx K - 2rx B + 3x s )] 

xHg\Ab B ,C 8 b B ,C 8 b K )F 8 (t 8b ) 
+ (1 - r)(l - r 2 )[(l + r)x B - rx K ] 

xHi c \y/\A»\b K , D 8 b B , D 8 b K )F(t 8c ) 
-[(1 - r 2 )((l - r 2 )(2 + x K ) ~ (1 - 3r)rr s ) 
-2r 2 x^(a; ft : - x B )] 



(1 - r 2 ){x B - x x )Ml , Bl = (1 - r 2 + x B )Ml , 

(1 - a^M 2 , D 8 2 = (1 - r 2 )x B M 2 , £ 8 2 = (1 - r 2 )a^M 2 , 

m&x(A,B 8 , l/b B , l/b K ) , 

max(A, C 8 , , 

max(y/\A!*\,D % ,l/bB,l/bK) , 
m a x(^\A^\,E 8 ,l/b B ,l/b K ) . 

F 8 , 

F 8 (t) = a s (t)c 8 (t) exp[-S(x B ,x K ,t, b B ,b K )} . 

>ns 

Hi a) (Ab K , B 8 b B , B 8 b K ) = Hf\Ab K , B 8 b B , B 8 b K ) , 

Hi b) (Ab B , C 8 b B , C 8 b K ) = K (Ab B )h'(C 8 b B , C 8 b K ) 

H { 8 c \^\A^\b K ,D 8 b B ,D 8 b K ) = 

K (J\A^\b K )h(D 8 b B , D 8 b K ) , A' 2 > , 
i|if? ) (V^&A:)MA»&B,A»&*r) , A' 2 < , 

Hi d \J\A^\b B ,E 8 b B ,E 8 b K ) = 

K (J\A^\b B )ti(E 8 b B , E 8 b K ) , A' 2 > , 
i|//,S 1) (^4^6 B )/ i / (£; 8 6 B , £ 8 &*) , A' 2 < 





are derived from Figs. 3(a)-(d), respectively. It is obvious that the above factorization 
formulas bear the features of Eq. fljjD . 

The exponentials exp(— S) appearing in M\ are the complete Sudakov form factors with 
the exponent 



S = s(x B P B , M + 2 / — 7 9 (a s (/x)) 

Jl/b B II 



/b B H 



+s(x K P£, b K ) + s((l - x K )P£, b K ) + 2 f ^ 7? (a s (/i)). 



The wave functions <pB and <pK* satisfy the normalization, 



4>i(x)dx 



I) 



2^6' 



(33) 



(34) 



with the decay constant fa, i — B and K*. The wave function for the K* meson with 
transverse polarizations has been derived using QCD sum rules [12], which is given by 

</> K . = %^(1 - i 2 ) * [0.267(1 - e? + 0.017 + 0.21£ 3 + 0.07£] , (35) 

with ^ = 1 — 2x. As to the B meson wave functions, we employ two models |2T| , [22 

N B x(l - x) 2 



M B + C B (l-x) 



4> { b\ x ) = n 'b^ x ( 1 ~ x ) ex P ( 



1 fx M B 



(36) 
(37) 



where Nb and N' B are normalization constants; while Cb and uj are shape parameters. 



IV. RESULTS AND DISCUSSIONS 

In the evaluation of the various form factors and amplitudes, we adopt Gf = 1.16639 x 
10 -5 GeV -2 , the flavor number n/ = 5, the decay constants fs = 200 MeV and fx* = 220 
MeV, the CKM matrix elements \V t * s V tb \ = 0.04, the masses M c = 1.5 GeV, M B = 5.28 GeV 
and M K * = 0.892 GeV, the B° meson lifetime r B o = 1.53 ps, and the QCD scale A = 0.2 



GeV [ID]. We find that, no matter what value of the shape parameter Cb is chosen, the 
model <j)]p in Eq. (|36| ) with a flat profile leads to results smaller than the CLEO data which 
gives B(B — > K*^) = (4.2 ± 0.8 ± 0.6) x 10 -5 0. In fact, the maximal prediction from (f) B , 
corrsponding to the shape parameter Cb = —M B , is about 3.0 x 10~ 5 , close to the lower 
bound of the data. 

On the other hand, using model (j)^ in Eq. (|37|) with a sharp peak at small x, we obtain a 
prediction much closer to the experimental data. It is indeed found that, as uj = 0.795 GeV, 
a prediction 4.204 x 10~ 5 for the branching ratio is reached, which is equal to the central 
value of the experimental data. If varying the shape parameter to both uj = 0.79 GeV and 
uj = 0.80 GeV, we obtain the branching ratios 4.25 x 10~ 5 and 4.14 x 10~ 5 , respectively. It 
indicates that the allowed range for u is wide due to the yet large uncertainties of the data. 
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The detailed contribution from each amplitude Mf is listed in Table II in the unit of 10~ 6 
GeV -2 . It is clear that Mf and Mf together give dominant contributions to the decay width. 
One also sees that Mf and Mf are smaller by an order of magnitude, while the amplitudes 
associated with 0$ are highly suppressed. In table II, it is interesting to note that Mf adds 
constructively to Mf, i.e., |Mf + Mf | 2 > |Mf | 2 . On the contrary, Mf is destructive to 
Mf. Due to this cancellating effect, the inclusion of O2 contribution only enhances the total 
rate by 2% (This result is basically consistent with the estimation obtained in JT3]]). This 
result is not sensitive to the choice of uj. We obtain the same enhancement for the total rate 
with uj chosen to be 0.79 GeV and 0.80 GeV, respectively. 

By fitting our prediction for the branching ratio BiB — > K*j) to the CLEO data, we 
determine the B meson wave function, 

<f> B {x) = 0.740079^(1 - x) exp \-\ ( - - VI , (38) 



2 V0.795 GeV 



which possesses a sharp peak at the low momentum fraction x. We stress that, however, 
Eq. (|38|) is not conclusive because of the large allowed range of the shape parameter uj. A 
more precise B meson wave function can be obtained by considering a global fit to the data 
of various decay modes, including B — > D^n(p) f23]| . Once the B meson wave fucntion is 



fixed, we shall employ it in the evaluation of the nonleptonic charmless decays. 

Finally, as mentioned in the Introduction, the authors of Ref. [|14|] found that the diagrams 
without hard-gluon exchanges dominate over Figs. 2(a) and 2(b) we have evaluated (only 
the operator O7 was considered in the calculation of the branching ratio B(B — > K*^) 
in |14[]) : the latter contribute at most 23% to the branching ratio, or 12% to the decay 
amplitude. However, in that analysis, a s is set to 0.2, which is even smaller than a s (M#) = 
0.23 evaluated at the B meson mass. We argue that such a small coupling constant is 
inappropriate, since the momentum flow involved in the decay process is most likely less 
than the b quark mass M&, say, roughly 1 ~ 2 GeV, which corresponds to a s ~ 0.4. The K* 
meson mass was neglected in [[14]], such that Fig. 2(a) does not contribute. In our analysis 
we did not make this approximation, and observed that the contribution from Fig. 2(a) is 
about Mk*/Mb ~ 1/5 of that from Fig. 2(b). Moreover, a flat B meson wave function 
corresponding to the shape parameter uj w 1.3 GeV was adopted in the computation of the 
higher-order contributions in [14]], which is far beyond uj ~ 0.4 GeV specified in ||22|| . If a 



sharper B meson wave function is employed, these contributions will be enhanced at least by 



a factor of 3. Note that the leading-order contributions considered in [ll] are less sensitive to 
the variation of uj in the wave function. Adding up the above enhancements, the amplitude 
for 0(a s )-corrections becomes approximately equal to the leading-order contribution. In 



this sense the analysis in |14[ does not seem to be self-consistent. 

In our approach, the momentum of the spectator quark in the B meson is parametrized 
into the minus direction, while the momentum of the spectator quark in the K* meson is 
parametrized in the plus direction. This parametrization is appropriate because of the hard 
gluon exchange. Note that the factorization formulas presented in Sect. Ill are constructed 
based on the diagrams with at least one hard gluon exchange. For example, the infrared 
divergences from self-energy corrections to the spectator quark are factorized into the B 
meson wave function, if they occur before the hard gluon exchange; and into the K* meson 
wave function, if they occur after the hard gluon exchange. Without hard gluons to dis- 
tinguish the initial and final states, factorization of self-energy corrections to the spectator 
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quark is ambiguous. Therefore, the diagrams without hard gluons do not appear in the 
regime of PQCD factorization theorems, and those with one hard gluon are indeed leading. 
If ignoring the validity of the factorization, we may proceed with evaluating the contribution 
of Oj from the diagram without hard gluons in the PQCD framework. It is trivial to obtain 
the hard subamplitude in momentum space, 

(1+ 0(1 ->*) ,. 



2M, 



B 



where the (^-function requires that the longitudinal momenta of the spectrator quarks in 
the B and K* mesons are in the same direction. Fourier transforming the above expression 
into b space, and convoluting it with the wave funcitons in Eqs. (35) and ( |3"8| ) and with the 
Sudakov factor, we derive the amplitudes 

M 7 5(0) = -M 7 P(0) = r (0) f 1 dx [ 1/A bdb<p B {x)<f) K ,{x) 



JO 

Without hard gluons, the momentum fraction x and the transverse extent b are equal for the 
B and K* mesons. We have set the hard scale t to 1/6 due to the lack of gluon momentum 
transfer. A simple numerical work on Eq. (EQ) gives M? {0) /TW = -2.51 x 1(T 6 GeV- 2 , 
which is of the same order as the contributions from 2 and Og. One of the reasons for the 



smallness of Mj , compared to the values obtained in |14[ , is the additional strong Sudakov 
suppression. Referred to Mf listed in Table II, is negligible. If including and 

M^°\ the branching ratio B(B — > K*^y) will increase by only 1.7%. 



V. CONCLUSION 

In this paper we have extended the PQCD three-scale factorization theorem to the 
penguin-induced radiative decay B — > K*j, which takes into account the Sudakov resum- 
mation for large logarithmic corrections to this process. We have included the non-spectator 
contribution from the current- current operator O2 besides the standard contribution given 
by magnetic-penguin operators 7 and 8 . It turns out that the contribution by 2 is neg- 
ligible due to certain cancellations. The contributions from 0§ and other operators in the 
effective Hamiltonian are also quite small. Finally, we have determined the B meson wave 
function from the best fit to the experimental data of BiB — > K* r y) ) which will be employed 
to make predictions of other B meson decays. 
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FIGURES 



FIG. 1. Contributions to the B — ► K*^ decay from the current-current operator O2. The 
diagram with a photon emitted from the other side of the charm-quark loop is not shown. 

FIG. 2. Contributions to the B — ► IT* 7 decay from the magnetic-penguin operator O7. 

FIG. 3. Contributions to the £> — > K*j decay from the chromo-magnetic-penguin operator Og. 

FIG. 4. Contributions to the i? — > K*7 decay from the strong-penguin operators O3, 04, • ■ ■ 
The dark square denotes insertions of the operators O3, 04, ■ ■ ■ C>6- 

FIG. 5. Contributions to the B — > K*7 decay from an C>2-insertion and a bremsstrahlung photon. 
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Table I. Hard subamplitudes obtained from Figs. 1-3. The quantities and are 
integrands of In and J 2 o — ho respectively, where the general integral I a b is defined in Eq. 
(10). 



Diagram H s 

4 (1 - r)(l - r 2 )x K [(l - r 2 + 2rx K + 2x B )A 4 + {rx K + 3x B )A 5 ] 



2 
7 (a) 
7 (b) 
8 (a) 
Os(b) 
Os(c) 
8 (d) 



3 x K x B m B + (k KT - k BT ) 2 

2r(l -r 2 ) 

[x K x B m 2 B + (k KT - k BT ) 2 ][(x B - r)m\ + k 2 BT ] 

2(1 -r 2 )[l + r + (1 -2r)x K ] 
\x K x B m 2 B + (k KT - k BT ) 2 ](x K m 2 B + k 2 KT ) 

(1 - r 2 + x B )(rx K + x B ) 
3[x K x B m 2 B + (k KT - k BT ) 2 ][(l - r 2 + x B )m 2 B + k 2 BT ] 

(2 - 3r)x K - x B + r(l - x K )(rx K - 2rx B + 3x B ) 
3[x K x B m 2 B + (k KT - k BT ) 2 ][(x K - l)m| + k 2 KT ] 

(1 + r)(l — t 2 )[(1 + r)x B — txk] 
3[(1 - r 2 )(x i? - x K )m% + (k^ T - k BT ) 2 ][(l - r 2 )x B m 2 B + k 2 BT ] 

(1 - r 2 )[(l - r 2 )(2 - x K ) + (1 + 3r)(2x y - x B )} + 2r 2 x^(^ - x B ) 
3[(1 - r 2 )(a; B - x K )M 2 + (k XT - k BT ) 2 ][(r 2 - l)x K m% + k\ T \ 
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Diagram H p 

4 (1 - r 2 )x K [{l - r)(l - r 2 ) + 2r 2 x^ + 2x B )A A + (r(l + r)x K + (3 - rx B ))i 5 ] 
Co — 



7 (a) 
7 (b) 
8 (a) 
Os(b) 
Os(c) 
Os(d) 



1 3 x K x B m 2 B + (k KT - k BT y 

2r(l -r 2 ) 



[x K x B m 2 B + (k^ T - k j B T ) 2 ][(x B - r)m 2 B + k BT ] 

2(1 - r 2 )[l + r + (1 - 2r)x K ) 
[x K x B m B + (k^ T - k BT ) 2 ](x K m B + fc|- T ) 

(1 - r 2 + x B )(rx K + x B ) 
3[x K x B m 2 B + (k KT - k BT ) 2 ][(l - r 2 + x B )m 2 B + k 2 BT ] 

(2 - 3r)x^ - x B - r(l - x K )(rx K - 2rx B + 3o;b) 
3[x K x B m 2 B + (k K r - k B T) 2 ][(^ - l)m 2 B + fc^ T ] 

(1 — r)(l — t 2 )[(1 + r)x B — txk] 
3[(1 - r 2 )(x B - x K )m 2 B + (k KT - k BT ) 2 ][(l - r 2 )x B m\ + k 2 BT \ 

(1 _ r 2)[ (1 _ r 2 )(2 + ^ _ (1 _ 3r)a . B ] _ 2r 2 rr g (a; Ar - Eg) 
"3[(1 - r 2 )(x B - x K )M 2 + (k KT - k BT ) 2 }[(r 2 - \)x K m\ + k 2 KT ] 



16 



Table II. The amplitudes M\ in the unit of 1(T 6 GeV~ 2 . 



M 2 A '/r 


Mf/r 


Mg/r 


-2.46 - u.m 


-140.14 - 143.31i 


-0.58 + 1.10i 


Mf/r 


Mf/r 


Mi'/To 


-1.21 - 11.05* 


140.14 + 143.3H 


0.54- 1.33i 
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